GLEON Kelly model test - first look
Dataset summary stats:
We have 62 lakes total. The table below shows the number of daily observations and date range for each lake. Lakes are arranged by number of daily metabolism observations.
| lakeName | n_daily_obs | min_day | max_day |
|---|---|---|---|
| PrairiePothole | 18 | 2019-07-05 | 2019-10-11 |
| PrairieLake | 19 | 2018-09-28 | 2018-10-16 |
| LittleRock | 20 | 2018-08-01 | 2018-10-02 |
| Balaton | 39 | 2008-06-13 | 2008-10-11 |
| Morris | 72 | 2013-05-21 | 2013-07-31 |
| Barco | 76 | 2019-04-12 | 2019-12-16 |
| Rotoiti | 82 | 2008-08-25 | 2009-07-23 |
| P1 | 92 | 2019-05-25 | 2019-08-24 |
| P8 | 92 | 2019-05-25 | 2019-08-24 |
| EastLong | 93 | 2013-06-13 | 2013-09-13 |
| Northgate | 93 | 2017-05-30 | 2017-09-14 |
| Cranberry | 95 | 2017-06-04 | 2017-09-14 |
| SkyPond | 96 | 2017-05-30 | 2017-09-25 |
| Ward | 99 | 2016-06-01 | 2016-09-11 |
| Bay | 100 | 2016-05-30 | 2016-09-10 |
| TheLoch | 100 | 2016-07-20 | 2016-10-27 |
| Bolger | 103 | 2017-05-31 | 2017-09-14 |
| Brown | 106 | 2017-05-30 | 2017-09-14 |
| Hummingbird | 108 | 2017-05-30 | 2017-09-14 |
| WestLong | 108 | 2017-05-30 | 2017-09-14 |
| Croche | 109 | 2016-05-10 | 2016-10-22 |
| Castle | 111 | 2019-06-12 | 2019-10-24 |
| Nastjarn | 113 | 2013-05-15 | 2013-10-13 |
| Crampton | 114 | 2013-05-21 | 2013-09-12 |
| Almberga | 120 | 2019-05-29 | 2019-09-25 |
| Trout | 131 | 2013-06-09 | 2013-11-18 |
| Lillinonah | 137 | 2014-05-06 | 2014-10-16 |
| Lillsjoliden | 137 | 2012-06-11 | 2012-10-25 |
| Mirror | 138 | 2003-03-20 | 2003-10-15 |
| Ovre | 141 | 2012-06-05 | 2012-10-23 |
| Onondaga | 142 | 2001-05-05 | 2001-10-27 |
| Erken | 145 | 2018-05-08 | 2018-11-06 |
| CrystalBog | 147 | 2008-05-30 | 2008-11-04 |
| Oneida | 147 | 2019-05-30 | 2019-10-27 |
| Struptjarn | 147 | 2012-05-31 | 2012-10-24 |
| Langtjern | 151 | 2013-06-12 | 2013-11-24 |
| Suggs | 152 | 2019-03-31 | 2019-12-12 |
| TroutBog | 153 | 2008-05-30 | 2008-11-10 |
| Acton | 156 | 2011-04-21 | 2011-10-02 |
| Simoncouche | 156 | 2016-05-19 | 2016-10-24 |
| Mangstrettjarn | 158 | 2013-05-10 | 2013-10-14 |
| Hampenso | 165 | 2007-04-19 | 2007-09-30 |
| Jordan | 171 | 2018-05-09 | 2018-10-31 |
| Sunapee | 175 | 2008-05-01 | 2008-10-30 |
| Utah | 182 | 2017-05-04 | 2017-11-19 |
| Taupo | 198 | 2017-10-01 | 2018-05-01 |
| Vedstedso | 201 | 2008-05-14 | 2008-11-30 |
| Schulzensee | 203 | 2010-05-07 | 2010-12-19 |
| FredriksburgSlotso | 211 | 2006-04-11 | 2006-11-21 |
| Harp | 212 | 2013-05-01 | 2013-11-30 |
| Gollinsee | 221 | 2010-05-05 | 2010-12-30 |
| StGribso | 226 | 2006-04-12 | 2006-11-23 |
| Pontchartrain | 238 | 2008-03-21 | 2008-12-31 |
| Taihu | 244 | 2007-10-09 | 2008-09-26 |
| Mendota | 246 | 2016-03-30 | 2016-12-01 |
| Mueggelsee | 292 | 2015-01-02 | 2015-12-29 |
| Rotorua | 310 | 2007-07-14 | 2008-07-09 |
| Vortsjarv | 319 | 2013-01-14 | 2013-11-30 |
| Feeagh | 344 | 2013-01-01 | 2013-12-31 |
| Annie | 356 | 2008-03-01 | 2009-02-27 |
| Kentucky | 357 | 2008-01-01 | 2008-12-30 |
| YunYang | 365 | 2016-12-31 | 2017-12-30 |
Distribution of log10 GPP
Distribution of log10 ER
Median vs. Mean
Time series of each lake
Plotted to give you a sense how the seasonality in GPP in some of these systems as well as what times of year we have data for each lake. 
Pull in lake metadata & join
Code not printed.
Maps & tables of global distrib.
Southern hemisphere and lakes in Asia not printed on a map, but note that they do exist!
| China | New Zealand | Taiwan | SUM |
|---|---|---|---|
| 1 | 3 | 1 | 5 |
Load calc possible?
Showing you how many lakes (n, as well as percentage) where we can get “true” load estimates from stream discharge and nutrient grab samples 
LSA distribution
Here I’m showing a few figures to give you a sense of the lake size distributions, both categorical and continuous. 
Hydrologic residence time distibution
Here I’m showing a few figures to give you a sense of the HRT distributions. These HRT estimates were either given to me by the folks sharing data, or I extracted them from the literature. Still missing a few lakes.
Pull individuals lake nutrients
Data wrangling, not shown.
Join nutrients to metab
Data wrangling, not shown.
In-lake nutrients
DOC, TP, TN distributions
We also have some other nutrient data for some lakes, but it is much less consistent. For instance, we have some absorbance data, but not always measured at 440nm. Some folks also supplied dissolved nutrients, and I believe all the particulate nutrients come from UNDERC lakes.
| name | Total number of lakes |
|---|---|
| abs250 | 6 |
| abs254 | 3 |
| abs280 | 6 |
| abs440 | 22 |
| chla_ugL | 23 |
| cond_uScm | 8 |
| DOC_mgL | 57 |
| DRP | 2 |
| NH4_ugL_N | 7 |
| NO3_ugL | 23 |
| NO3_ugL_N | 5 |
| orthoP_ugL | 6 |
| PO4_ugL | 7 |
| POC_mgL | 8 |
| PON_mgL | 8 |
| PP_ugL | 9 |
| secchi_m | 4 |
| SRP_ugL | 17 |
| TDN_ugL | 8 |
| TDP_ugL | 9 |
| TN_ugL | 53 |
| TP_ugL | 55 |
GPP vs. in-lake nutrients
In lake DOC vs. TP
GPP vs. lake size
Lake DOC vs. lake size
zMix vs. lake DOC
Modelled zMix vs. calc zMix
At our meeting last week someone mentioned something about comparing the Kelly model calculated zMix values (Equation 7 from the paper) with modeled zMix values. This is how it looks. We don’t have a lot of chla or POC data, which might be one reason this looks off. 